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Abstract 

We consider nonequilibrium transport in a simple chain of identical mechanical 
cells in which particles move around. In each cell, there is a rotating disc, with 
which these particles interact, and this is the only interaction in the model. It was 
shown in [1] that when the cells are weakly coupled, to a good approximation, 
the jump rates of particles and the energy-exchange rates from cell to cell follow 
linear profiles. Here, we refine that study by analyzing higher-order effects which 
are induced by the presence of external gradients for situations in which memory 
effects, typical of Hamiltonian dynamics, cannot be neglected. For the steady 
state we propose a set of balance equations for the particle number and energy in 
terms of the reflection probabilities of the cell and solve it phenomenologically. 
Using this approximate theory we explain how these asymmetries affect various 
aspects of heat and particle transport in systems of the general type described 
above and obtain in the infinite volume limit the deviation from the theory in [1] 
to first-order. We verify our assumptions with extensive numerical simulations. 

1 Introduction 

In this paper, we study nonequilibrium transport in a class of 1-D Hamiltonian 
systems consisting of free noninteracting point particles of mass m moving inside 
a chain of identical cells, which are like "chaotic billiards". Each cell contains 
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a device (called a tank in [1]) that interacts with the particles by exchanging en- 
ergy with them. Even though the particles do not interact among themselves, an 
effective interaction is mediated by the tanks. As pointed out in previous stud- 
ies [2, 3], this mediated interaction among the particles allows such models to 
reach thermalization. These models are therefore genuine many-particle interact- 
ing Hamiltonian systems [4]. Moreover, the details of the interaction between 
tanks and particles are not important. However, since our aim is to study model 
systems with realistic microscopic dynamics we shall ask that this interaction sat- 
isfies some general conditions: The system has to be time-reversible and Hamil- 
tonian, although conservation of phase space volume is seemingly enough. 

An earlier study, [1], went some ways in explaining the energy and particle 
profiles in terms of a stochastic approximation to a Hamiltonian model. This class 
of models was derived, in turn, from work [2, 3] in which a Lorentz gas with 
rotating discs was considered. 

The stochastic approximation used in [1] assumed that in their evolution, the 
particles exit on both sides of each cell with equal probability at some fixed rate. 
Out of equilibrium, the exit rates will change from cell to cell and this leads to the 
effective transport of particles and heat. It was found that, to lowest order, due to 
the gradient character [5] of the system, the profiles for the rates at which particles 
and energy are transported among neighboring cells interpolate linearly between 
the values imposed by the baths at the ends. From this, energy and particle density 
along the chain were computed explicitly. They are generally not linear. 

The equal probability to jump to the left or to the right leads to a simple ran- 
dom walk process, without any memory effects. This assumption is strictly valid 
only when the cells are weakly coupled, i.e., when the size of the openings con- 
necting neighbouring cells is very small. However, away from this limit case, the 
memory effects cannot in general be neglected and the symmetry assumption is 
likely to fail. 

In this paper, we study the consequences of including these dynamical mem- 
ory effects. The price we pay for doing this is the lack of a simple stochastic for- 
mulation to describe the energy exchange from cell to cell. Instead, at the steady 
state we propose a set of balance equations for the particle number and energy in 
each cell in terms of the reflection and transmission probabilities for particles and 
energy. We specify a phenomenological law for these probabilities that allows us 
to obtain an approximate expression for the steady state. 

Our findings for the class of 1-D Hamiltonian chains that we consider can be 
summarized as follows: 
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1 . We elaborate a phenomenological theory for the reflection and transmission 
probabilities for particles and energy. 

2. Close to equilibrium, one can model the Hamiltonian system by persistent 
random walks. In this case the system remains of gradient type and thus, 
the theory of [1] is still applicable with corrections of the order of with 

being the size of the system. 

3. Far from equilibrium, we demonstrate that there is indeed a non-negligible 
dependence on the local gradients, and the corresponding profiles can be 
shown to obey a non-linear differential equation that goes beyond linear 
response regime. 

For our systems this asymmetry can be seen as a result of two different con- 
tributions: one which is purely geometrical and another which has a dynamical 
origin. For any scattering billiard (like the model cells) the reflection probability 
will in general be different from the transmission purely due to the geometry of 
the cell. When some source of interaction among the particles is considered, these 
coefficients will depend in addition on the local thermodynamical fields (namely, 
density and temperature) of the particles in the cell. When, in addition, the cell 
is subjected to an external thermodynamical gradient, the reflection and transmis- 
sion probabilities will be different at the left and at the right of the cell. In this 
situation the dynamical asymmetry is emphasized. 

Close to equilibrium, the geometric component dominates the dynamics, and 
then the asymmetry is independent of the local thermodynamical fields and thus, 
uniform along the chain. Furthermore, in this simpler case, the bias of the reflec- 
tion/transmission coefficients is the same on both sides of any cell. The corre- 
sponding stochastic model is a persistent random walk process like that studied 
in Lorentz gases [6]. After introducing the model in Section 2, in Section 3 we 
sketch the stochastic approximation on which the results of the subsequent sec- 
tions are inspired on. In Section 4, we study the case of constant asymmetry and 
show that one recovers the gradient property observed in [1], albeit with a slight 
correction in the slope of the jump rate profiles. 

In the final sections, we analyze how the presence of non uniform local fields 
changes this picture. Clearly, since the tanks induce an interaction of the particles, 
the reflection (or transmission) probability depends on the local thermodynamical 
fields. We will distinguish two contributions to the dynamical asymmetry: one 
that depends on the mean values of the fields and another one that depends on the 
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local gradients of the fields. These two effects are added to the purely geometrical 

contribution of the corresponding non-interacting system. 

In order to judge the importance of the asymmetry phenomenon, we take fixed 
external fields and let the number of cells be sufficiently large. This is the stan- 
dard limit taken in studies of the Fourier law [7, 8]. In our setting, the system 
is no longer of gradient type, and we obtain expressions for the transport equa- 
tions which go beyond the standard linear regime. While these equations remain 
approximate, we verify numerically their validity even very far from equilibrium. 

2 The Mechanical Model 

We consider a gas of noninteracting point particles of mass m that move freely in- 
side a one-dimensional chain composed of identical two-dimensional cells, ar- 
ranged horizontally. Each cell is connected to its left and right neighbors through 
two openings of size 7. 

Inside each cell there is a mechanical "tank" capable of storing a certain 
amount of energy. When a particle hits the tank, some energy exchange takes 
place, and this is the only interaction in the model. 

While the details of this interaction are largely irrelevant for the discussion in 
this paper, our simulations have been done for the following precise setup which 
we call the rotating disc model (RDM). Each cell has the geometry described in 
Fig. 1. The energy tank is modeled by a freely rotating disc which is pinned at the 
center of the cell. Its energy is only rotational. When a particle hits the disc, it 
exchanges its tangential velocity v^. with the angular velocity u of the disc, while 
the normal component is reflected [2]: 

i^'^Vt, v't^u, v'^^-v^. (2.1) 

Note that while the particles do not see each other, they effectively interact through 
their collisions with the tanks. Thus local equilibrium can be reached [9]. The 
model is by no means free. 

The chain is connected at both ends to two reservoirs of particles through 
openings of the same size 7. The reservoirs are idealized as infinite chambers 
containing an ideal gas at a certain density n and temperature T. When we study 
non-equilibrium effects, n and T will differ at both ends. 

Another variant of this model has the same fixed walls, but the disc is replaced 
by a rotating needle, with elastic collision, conserving total energy (and total an- 
gular momentum). 
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Figure 1: Upper left: The geometry of a single cell for the rotating disc model. The 
boundary of the cell is made of four arcs of circles of radius R^. The radius R of the 
rotating disc is chosen so that no trajectory can cross the cell without undergoing any 
collision. The arrows represent the incoming (j) and outgoing (J) rates of particles that 
cross the exits of width 7. The rates q and Q for the energy are not shown. Lower right: 
A chain of 5 cells connected to two baths. 



These models are among a class of models satisfying the following minimal 
assumptions: 

• Time reversibility: The trajectory of the system in phase space across a 
scattering event must be time reversible. While this condition is not re- 
ally necessary for our derivations, our aim is to model realistic macroscopic 
situations for which microscopic time-reversibility is believed to hold. 

• Conservation of phase-space volume: This condition is not as strong as 
requiring the dynamics to be Hamiltonian, and seems to lead to similar re- 
sults. In particular, the rotating disc model satisfies this condition but is not 
Hamiltonian: Its collision rule preserves phase space volume but is not a 
canonical transformation. 
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We call the models satisfying these two assumptions mechanical. 

We end the description of the model by a more detailed account of the heat and 
particle baths. From the left reservoir at density and temperature T>, particles 
are injected into the system at a rate Observe that the mean energy of the 
particles injected into the system is not T> but |T> (see e.g., [1]). Thus, the left 
reservoir injects energy into the system at a rate q-^ given by 

3 . ^ 

Analogously, the right reservoir injects particles and energy at rates j< and g< 
respectively. Each injected particle will be eventually re-absorbed by one of the 
reservoirs and this happens when it crosses into that reservoir. All information 
about that particle is then discarded. 

As the reservoirs consist of an ideal gas, the injection rate j is a function of 
the density n and temperature of the reservoir given by 

j oc 7nrV2 _ (2.2) 

Furthermore, the chemical potential at the reservoirs can be written in terms of the 
injection rate as 

^'-Tlog(J0,^ (2.3) 

with Aq some constant. 

We end the section with some comments on the choice of the model: Our 
choice of the reflecting boundaries of Fig. 1 is to guarantee the desirable ergodic 
and mixing properties of chaotic billiards. However, we have not succeeded in 
showing rigorously that the cell is either ergodic or mixing. While the geometry 
of the rotating needles model (RNM) is similar to that of the RDM, there are 
two main differences: First, the RNM is not only mechanical but Hamiltonian. 
Second, in the RNM model, a particle can in principle cross many cells without 
hitting any needles or boundaries. This may produce logarithmic corrections to 
transport, which are beyond our study. 

3 An Approximate Stochastic Model 

In [1], a stochastic approximation of the mechanical model was considered. In 
this approximation it was assumed that what happens in one cell is independent 
of the state of the other cells. The particles perform a random walk among the 
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cells, mixing their energies with the tanks according to rules which mimic the 

mechanical system described above. Once the particle has mixed its energy, it 
waits (until an exponential clock rings) and then jumps with equal probability 
to the left or to the right. With these assumptions it was found that the rates at 
which particles and energy jump among neighboring cells are linear functions of 
the position and from this, energy and particle density along the chain can be 
computed explicitly under very general assumptions. However, note that these 
assumptions are valid when the cells are weakly coupled and dynamical memory 
effects appear when this is not the case. 

The purpose of this paper is to study the consequences of these memory effects 
on the transport properties of the mechanical models that we consider here. The 
price we pay is that the stochastic model, used in [1] is not longer appropriate. We 
briefly sketch the essentials of this stochastic approximation that will serve us as 
a guideline to our analysis (we refer the interested reader to [1] for more details 
about this). 

The stochastic description is based on particles, carrying energies, and the 
discs (as in [1]). Each state of the cell is represented by a point in Q = U^r^™^ 
where = {(x^, . . . , a;„J, (p^^, . . . , p„^), y} and m is the number of particles in 
the cell at some given time, are their energies, G { — 1,1} specifies the side 
on which the particle enter the cell and y is the energy stored by the disc. For a 
chain with N cells, the phase space is Q^^^ — QxQx • • • xQ with N factors. The 
reservoirs at the ends absorb particles, and eject them at a rate resp. (for left 
and right) and with an exponential distribution of energy (with mean temperatures 
T> and T<). 

The stochastic process is as follows (see [1] for more details): When a particle 
entered a cell k from the left, it will wait for an exponentially distributed time 
(whose mean may depend on its energy), it mixes energy with the disc, and leaves 
the cell (after another exponential waiting time) with probability k 
and with probability 1 — 0;^ it to the right. 

This stochastic description leads, in the steady state, to the balance equations 
for rates of particle injection (see (3.1) below). 

However, while this description is adequate for particle flux, it is too detailed 
for the energy flux. Instead of trying to formulate a stochastic process for both 
particle and energy fluxes (see [11] for an example), we elaborate a phenomeno- 
logical theory for the physical parameters involved. We show that this information 
is enough to obtain precise results on the average properties of the steady state of 
this type of systems. 

It is convenient to distinguish the local incoming rates from the outgoing rates. 
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We denote by Jj^ ^ and Qi^ ^ the rates at which particles and energy exit the k- 

th cell to the left and by ^, ^, those to the right. The exit rates are then 
simply determined by the corresponding incoming rates Jl/r.^ Qh/R.k^ by 
reflection a^i^, ^ coefficients which will satisfy the balance equations (for the 
cell k): 

JL,k = Ct^fc iL,fc + (1 ~ Ck;R,fc) jR,k ) 
jR,k = (1 ~ CtL,fe) iL,fc + CtR,fc jR,k ) 

o o (3-1) 

QL,k = "L,fc 9L,fc + (1 - «R,fc) gR,fc , 
<3R,fc = (1 - a?,fe) 9L,fe + aR,fe QR,k ■ 



Note that even for (left-right) symmetric cells, out of equilibrium the reflection 
coefficients at the left and right of the cell are not necessarily equal. The balance 
equations account for the conservation of energy and particle number, namely 
QL,k + QR,k = (lL,k + (lR,k and JL,fc + JR^fc = jL^fe + jR^fc. Since an average the 
steady state is assumed, mean rates appear in (3. 1). 

We define the coefficients by empirical probabilities Jlr ... by (omitting the 
index k for better legibility) 



J _ Jll j _ Jrr 

Jl Jr 

Q 1 , ?LL ~ Q.LR Q 1 , QrR " QrL 

ttt =7^ H ^ , ttp =7^ + 



(3.2) 



2 2q^ ' ^ 2 2gR 

where, for example, Jlr is the rate of particles leaving on the right, which entered 
on the left, and ^lr is the mean energy carried out by particles which entered on 
the left. If we denote p(E'\E) the conditional probability that a particle entering 
with energy E leaves with energy E' and by au^(E', E) the probability that it 
leaves on the left when it entered on the left, then one can think of ^s 

gLL = j dE'dEE'a^^(E',E)p(E'\E), 

where the complicated dependence on the state of the system is neglected (that is, 
we do not consider the density and the other particles in this formula). Instead we 
will use a phenomenological description for the reflection coefficients. Perhaps it 
is useful to note that with the above conventions. 



Jll = J dE'dEa^^(E',E)p(E'\E). 
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The definition for the a is canonical, but for the we chose a more com- 
plicated expression: It simultaneously preserves total energy conservation, but 
allows for an energy change during scattering. The reader should note that the 
quantities a, J, Q are mean values averaged on the (fluctuating) steady state. An 
important aspect of the present work is to check that this approximation still cap- 
tures the essentials of transport of heat and particles. 

To further simplify the discussion, we make the approximate assumption that 
the a§ (for all X E {J, Q} and Y G {L, R}) only depend on the incoming fluxes. 
Due to the mechanical nature of the models considered here, the time scale is a 
free variable leading to the scaling relation 

«0l, Jr, Ql, Qr) = tt(AjL, ^Jr, A^^l, A^^r) , (3.3) 

for all A > 0. Therefore, a only depends on 3 ratios. It will be useful to distinguish 
the contribution to a that arises out of equilibrium from the contribution that only 
depends on the mean values of the fields. Accordingly we write 

and require that £(j^/^/g^/^, 0, 0) = 0, i.e., that e vanishes at equilibrium. The 
term Uq, which has a purely geometric origin, describes in turns those aspects 
which hold at thermal equilibrium and the term j^^'^/q^^'^, with q — (qi^ + ?r)/2 
and j = (jl + Jr)/2, is proportional to the density (see Eq. (2.2)). 

Note that ay = \ corresponds to the simple symmetric random walk consid- 
ered in [1]. The case when the a§ are independent of the local fields, but different 
from I corresponds to a persistent random walk and will be studied in the next 
section. More general and realistic laws will be investigated in the last sections. 

Also note that the details of the dynamics, in particular those that depend on 
the specific model, are encoded by Eq. 3.4. To show the validity of our approach 
we will proceed as follows: We first introduce the law (3.4). This closes the 
balance equations for the energy and particle number. We next determine the 
thermodynamical profiles and currents. Finally we compare these "theoretical" 
profiles to those obtained with numerical simulations of the mechanical model. 

4 Solution of the Stochastic Model Near Equilibrium 

In this section, we discuss profiles in the infinite volume limit with fixed boundary 
conditions, under the approximating assumption that every cell is at local equilib- 
rium. This means that we neglect the term e in Eq. (3.4) and thus, ckr j^. = ifc — 
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ckq and similarly for the a^. Furthermore, we make the assumption, that and 
are independent of the densities found in the chain. In this approximation the 
reflection coefficient is a constant independent of the thermodynamical fields and 
thus, independent of position. The exact range of validity of this assumption will 
be discussed in Sec. 5. Our approximations mean that we are dealing with the 
case of a persistent random walk, and the current section, while "well-known" in 
principle, serves as a first check of the validity of our approximations. Once this 
check has been done, we will study in more detail the corrections to the persistent 
case. 

Since the geometric contribution is the only one remaining in (3.4), we refer 
to this case as the geometric approximation. Note that this contribution will 

lead in general to a reflection coefficient different from |. This approximation 
corresponds to a persistent random walk in which the probability for the walker 
to move forward is different from the probability to move backwards. 

We discuss the total ejection rates = J^^ + J^^ and = Qj^j. + Qs_^k- 
The boundary conditions of the problem are simply 

j> = Jl,! , j< = jR,N, Q> = <1l,i , Q< = <1r,n ■ (4.1) 

From now on we will use rescaled variables ^ = k/(N + 1) E [0, 1], and let 
J(0 — Jfc and Q(^) — Qy.. For the sake of simplicity, we will only write and 

Lemma 4.1 Under the assumptions made above, the particle and energy ejection 
rates satisfy 



( 1 - 2a\l - 20 
1 + (A^ - l)a^ 



AO = 20> + e • Aj) + \ , ... T Aj , (4.2) 



where Aj = — j^, and 
with Aq — — q>. 

Remarks: The first term in the r.h.s. of Eqs. (4.2) and (4.3) corresponds to the 
result obtained in [1]. The last terms account for the correction when a 7^ |. 
These corrections disappear in the infinite volume limit. The integral over ^ of the 
correction vanishes, because of energy and particle conservation. Also, note that 



Solution of the Stochastic Model Near Equilibrium 



11 



as a is a constant, the profiles are linear (even when a ^ \), because the system 
is still of gradient type. 

Proof: The proof is just a calculation. We will concentrate on the case of particle 
rates. This calculation starts from (3.1). Using the identities Eq.(4.1), we can 
solve (3.1) for the remaining J and obtain fox k = 1, . . . ,N , 



(N -k + l)a . l + (k-2)a . 
~ 1 + (A^ - l)a'^>^ 1 + (N- l)a-^< ' 
_l + (N-k-l)a . ka . (4.4) 



1 + (A^-1)q; 1 + (A^-1)q;- 



Then, summing the R and L terms, we get 

_ 1 + 2(iV - k)a . 1 + 2(k l)a . 



Rearranging terms one immediately obtains (4.2). The case (4.3) is handled simi- 
larly. 

Therefore, if the value for a is known, the profile of the ejection rates can 
be obtained from the solution (4.4) for any nonequilibrium state, as long as a is 
independent of the external parameters. 

4.1 Thermodynamical profiles 

In order to identify temperature and particle density, we use the techniques of 
[1]. In that paper, it was observed that in equilibrium, a single cell has a Poisson 
distribution of particle number with mean n and a canonical distribution of the 
total energy in the cell (particles and disc) with a temperature T. Furthermore, 
(see Proposition 4.1 in [1]) one has universal relations between these quantities 
and the properties of the reservoir, namely 

T -19. 
3 J ' 

J3/2 
^ - V0-QTJ2 ■ 

The constant t]q depends only on the geometry of the cell and the size of the 
opening, but not on the shape 
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where the area is that of the cell minus the disc and where I7I is the size of the 
opening between adjacent cells. This identification is unique provided the system 

is ergodic' 

Assuming local equilibrium at any cell k in the system, the formulas above 
generalize immediately to predictions of the profiles for temperature and particle 
number (with ^ = k/(N + 1)) 

T(a = - ^ • (4 7) 

3 J(0 ' ^ ^ 

Note that while the jump rates Q and J have linear profiles in the current approx- 
imation the profiles of T and n are generally nonlinear. In Sect. 6 we will study 
in detail the deviations from linearity of Q and J. 

To test these results we performed out of (but close to) equilibrium simulations 
for a chain of 20 RDM cells with I7I = 0.08, = 1.15 and a disc radius of 
R = 0.0793 (ensuring that no particle can cross a cell without being scattered). 
The choice parameters used in the simulations are a compromise between good 
mixing within each cell, and the speed of convergence to the stationary state. 

We checked in our simulations for the Hamiltonian model that the profiles 
of the outgoing rates are described by (4.2) and (4.3). In Fig. 2 we show the 
temperature and density profiles for two different simulations, in which the baths 
were set to j> = 10, j< = 12, T> = 1000 and T< = 1100 (circles), and j> = 12, 
j< = 10, T> = 1000 and T< = 1300 (triangles). 

To compare the measured profiles for T and n with our analytical expressions 
we have first numerically computed the reflection coefficients and for each 
cell. For both simulations we found that both coefficients are constant (to within 
numerical accuracy) along the chain with a mean value of — 0.5585 ± 10~^ and 
a'^ — 0.5609 ± 10"^. We insert these values into (4.2)-(4.3) and obtain estimates 
for the profiles of the outgoing rates J(0 and Q(0- We then use these estimates in 
(4.7) and (4.8) to obtain the profiles for the energy T and the number of particles 
n, as predicted by our theory. These are shown as solid lines in Fig. 2. The fit with 
the directly measured profiles for T and n is excellent. Note that these profiles are 
not expected to be linear [1]. This non-linearity has already been verified against 
the theory in [1], but here we study rather the effect of the asymmetry. At this 



^The defimtion of t/q differs from that of [1] because of different normalizations of the kinetic 
energy. 
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Figure 2: Comparison of the measured profiles witii tiiose predicted by the Eqs. (4.7) 
and (4.8), using (4.2)-(4.3). The measured profiles were obtained as an average over 400 
different realizations of the RDM for 20 cells. The energy profile T{^) (a) and the particle 
density profile n{^) (b) are shown for two different simulations with = 10, = 12, 
r> = 1000 and T< = 1100 (circles) and j> = 12, j< = 10, T> = 1000 and r< = 1300 
(triangles). For all data the error bars are smaller than the symbol size. The solid lines in 
each panel correspond to the analytical profiles of Eqs. (4.7) and (4.8) for = 0.5585 
and = 0.5609 as explained in the text. 



order of approximation, that is, in the case of the persistent a's the difference 
between the two approaches is only of order and is not yet visible in the 
simulations. 

4.2 Macroscopic currents 

The heat and matter local macroscopic currents can be easily derived from the 
balance equations. The current of particles between the k-th and k + 1-st cells is 
defined as ip^ = J^^ — j^^ or, in terms of outgoing jump rates as ip^ = J^^ — 
JL,k+i- Using Eq. (4.4) we obtain 

1 - 

1 + (A^ — l)a^ 

where the result is of course independent of k. In an analogous way, = Qr ^ — 
Ql./c+I' ^iid we obtain for the heat current 

1 - 

m^ = -Aq. (4.10) 
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In these variables, the currents are linear functions of the external gradients. 

In the numerical simulations described above, we also computed the heat and 
matter currents obtaining an average value of: Lp^ = —0.077 ± 0.002 and ip^ = 
—185 ± 13 (for the experiment in circles in Fig. 2) and (p^ = 0.08 ± 0.008 and 
(p^ — —56 ± 4 (for the experiment in triangles of the same figure). If instead we 
insert the numerically computed values for and into Eqs. (4.9) and (4.10) 
we obtain (p^ = —0.076, v?^ = —180.8 (circles) and (p-^ = 0.076, = —56.5 in 
perfect agreement with the numerical experiment to within numerical accuracy. 

Therefore, we conclude that in cases in which the a probability can be taken 
to be constant along the chain, the conclusions of Lemma 4. 1 describe remarkably 
well the transport properties of the disc model. Similar results were found for the 
needle model, and we conjecture that this extends to many similar models. 

The validity of the linear transport equations (4.9) and (4.10) depends crucially 
on the property that q;l = aR. For high gradients, this property is violated as we 
show in Sec. 6. In that case, the Fourier law does not hold, but the flux still scales 
asl/N for large systems. The reason is that local gradients appear in the diffusion 
constants. 

We remark finally that a (and thus the details of the cell) determines the nature 
of the macroscopic transport: If a = the transport along the chain is ballistic 
as the currents do not scale with the size of the system N. On the other hand, if 
q; = 1 the chain behaves as an insulator. For any other value of a the transport is 
normal with well defined transport coefficients. 

5 Validity of the Geometric Approximation 

In the previous section, we assumed that the reflection coefficient a, while not 
necessarily equal to |, neither depends on the particle density inside a cell nor 
on the strain acting on it. In other words, «^ ^ = ^ = ctQ, this value being 
determined by the geometry of the cell only (and similarly for a^). 

In this section we discuss this assumption and determine the range of param- 
eters for which it applies with high precision. The discussion will also clarify 
the choice of the decomposition (3.4). The profiles obtained when and are 
different and vary along the chain are discussed in the next section. 

The reader should first note a simple fact: When the discs (tanks) in the cells 
are not allowed to rotate, and the reflection is specular, then the trajectories are 
independent of the energy and of the particle density. In that case a is clearly 
independent of all external parameters (and — a*^ since the particles are the 
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only energy carriers). Therefore, any variation of a has its origin in the effective 
interaction of the particles. 

The interaction changes the dynamics inside the cell, but does not always in- 
duce a variation of a. We first consider an ideal case: particles are always trapped 
in the cell for a very long time. They collide with the tank many times, and for- 
get the side from which they entered the cell, as well as the energy they had at 
that time. In this case, the thermalization of the particles is perfect due to the 
interaction and no asymmetry appears. 

The geometric approximation loses its validity due to two related memory 
effects: 

• The distribution of energy in the cell is not uniform, and thus thermalization 
is only approximate 

• Particles retain a memory of the energy they had when entering the cell, and 
this effect is energy-dependent. 

In order to test these properties, we first performed equilibrium simulations. 
These will isolate the geometric contributions to the a's. This allows us, in the 
next section, to identify those contributions which are typical of non-equilibrium, 
justifying thereby the decomposition of Eq. (3.4). In Fig. 3 we show the depen- 
dence of the reflection probabilities ckq and q;§ on the density of particles n for a 
single cell (for the RDM). Injection rates of the reservoirs are chosen according to 
(2.2) in order to obtain the desired densities. We observe the following: 

• For low densities (n < 0.3) and high densities (n > 10), the a is practically 
constant. Therefore, if a whole chain of cells is in either of these regimes, 
we can model it with a constant a. If in addition, the local gradients are 
small, then the system is (close to) gradient type, and is well described by 
the approximations of Sec. 4. 

• For intermediate densities, a increases linearly with the density. This is a 
new kind of regime (not of gradient type), which we discuss in more detail 
in the next section. 

As discussed in Sec. 3, the distinction in Eq. (3.4) between the geometric 
contribution ckq and the dynamical e, is that while e depends on the local gradi- 
ents of the thermodynamical field (and is zero at equilibrium), aQ will depend at 
most on the mean fields. Therefore, we can expect that a constant reflection co- 
efficient will always remain a good approximation close to equilibrium, namely 
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Figure 3: Dependence of the equilibrium reflection probabilities ajj (squares) and 
(circles) on the density of particles in the cell n for the RDM with 7 = 0.08. Note 
that aQ saturates to a constant value at low and high densities. At intermediate densities 
0.3 < n < 10, Oq grows linearly with density as corroborated by the linear fits (dashed 
lines). The fits are ag(n) = 0.56048 ± 2 • 10"^ + (1.29 • lO"'' ± 5 • 10-^)n, and a^(n) = 
0.55813 ± 2 • 10"^ + (8.71 • 10"^ ± 5 • 10-^)n. 



when AT/T ^ 1 and An/n ^ 1. In that case, the equations of the Sec. 4 hold, 
regardless of the average density of the system. 

The reader should observe that and «q differ, while for a system without 
memory effect (such as the Lorentz model) one would expect equality of these 
quantities. This difference is not an artifact of our simulations but a property of 
mechanical models of the type we consider. Looking at Fig. 1, and the collision 
rule (2.1), the reader will realize that fast particles have a higher probability to exit 
after only 1 collision with the disc than slower ones, because the angle of reflection 
tends to be smaller for the fast ones. Therefore the memory effect depends on the 
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Figure 4: Difference Aa = — as a function of AT/T for (squares) and 
OL^ (circles). Each symbol corresponds to a simulation for the single cell RDM with a 
geometry as in Fig. 2 but |7[ = 0.08. In all experiments the external gradients were 
fixed to j> = j< = 1, AT = T< - r> = 100 and \{T^ + T<) = T. The dashed 
lines correspond to Aa ~ {l^T /T)^/"^ . In the inset and at the left (plus) and right 
(crosses) boundaries are shown. Note that Aa ^ on approach to equilibrium, while for 
very large gradients reaches the Lorentz gas limit indicated by the hashed triangle (see 
the text for a discussion). 



individual energies and not the mean, and this accounts for the difference above. 
We have checked that by placing non-rotating discs between the openings and the 
turning disc, the effect is reduced. 

We now turn our attention to nonequilibrium effects, which are described by 
the functions Sy of Eq. (3.4). We performed simulations for the single cell RDM 
in which the external gradients were kept fixed but the mean temperature was 
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changed (T^^T- AT/2 and T< = T + AT/2). In the inset of Fig. 4, we show 
the left and right reflection coefficients, aj^ a^, as a function of AT/T. As 
the parameters approach equilibrium values, the difference |«l ~ <^rI ^^^^ to 
decay to zero (at a rate (AT/T)^/^). The value that al^ and take in equilibrium 
is the constant geometric asymmetry (which for these parameters are ckq ~ 0.558 
and ~ 0.560). At the other extreme, since the temperature cannot be nega- 
tive, AT/T < 2. The difference Aa is maximal at AT/T = 2. Interestingly, 
^[(AT/T = 2) corresponds to the reflection coefficient of the Lorentz gas, i.e., 
to the limit in which the central disc does not rotate. 



6 Far from Equilibrium 

In Section 4 we have assumed that the system is at local equilibrium, and that 
the particle density varies in a range where the are essentially constant. In 
this section, we discuss the case when both these assumptions are dropped. This 
means that we take into account the dependence of the a on the density and on 
the local gradients of the thermodynamical fields. 

Far from equilibrium, the distribution of the local fields will be different from 
cell to cell and thus, the probabilities will depend on the position along the chain 
a = q;(0- We recall the set of equations Eq. (3.4) for a one of the functions ay- 



Oiik, Jr, Qi^, Qr) = OtQ 



It is understood that ckq and e denote functions which are different for the various 
ay, but have the following common properties: By definition, e must vanish at 
equilibrium 

e^(a;,0,0) = . 
Furthermore, the left-right symmetry of the cell implies: 



y, z) = £^(x, -y, -z) . (6.1) 

We next discuss specific properties of e as they appear from numerical simu- 
lations for the RDM. It turns out that the e§ not only vanish on the submanifold 
(x, 0, 0), but, to a very good approximation they vanish for those points where 
the temperatures Tl and are equal. Since the temperatures are functions of 
the injection rates of particles and energy this means that e§ vanishes on the sub- 
manifold of {jL^iR) ?r} where Tl = T^. While we have no proof of this 
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Figure 5: Contour density plot of Ol ^ function of the local gradients in j and T. The 
obtained values correspond to the one-cell RDM with an opening of 7 = 0.08. Note that 
the contour line for Tl = Tj^ is horizontal, justifying Eq. (6.2). 



observation, it can be understood by observing that for = the tank has the 
same temperature as all the particles, and therefore we are in presence of purely 
geometric phenomena which are already captured by the function alone. 

We expand Cy to first order, and in view of the above information, a very good 
approximation is given by assuming that e has the form: 

^r + Tl Tr + Tl Jr + Jl 

where Sq and are constants depending on the choice X, Y . We have numerically 
corroborated Eq. (6.2) for the RDM. In Fig. 5 a contour density plot for a]^ is 
shown as a function of the local gradients. The linearity of the contours shows 
that is ^ linear function of the local gradient in j whose slope and intercept 
depends on the local gradient of temperature. The same behavior is found for the 
other coefficients. 

The reader should note that for fixed external parameters (for the heat and 
particle reservoirs) and for long chains, the local gradients per cell scale like 
and thus, one expects the e to be both proportional to 1 /N and to the force fields. 
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The question then is whether the effect of the asymmetry will disappear in the 
infinite volume limit or not. 

We next show how, in principle, the above approximations lead to a closed 
system of equations for the profiles, and we then check, for one particular case 
that the equations indeed describe the profiles found in numerical experiments, 
and account correctly for the deviation from [1]. 

Since the local gradients are naturally expressed in terms of the injection rates 
we use the identification between injection and ejection rates to write Eq. (3.1) as 

jR,k-l — 'A.,k jL,k + (1 ~ O^R./fc) jR,k 1 

jL,k+l = (1 - «[,ik) Juk + «R,fe jR,k > (6-3) 



with k — 1, . . . , N and boundary conditions j^i — j> and = j<. 

Using the definition of the particle current ip^ = j^k+i ~ Jr,*;' one can eliminate 
the (or the ji), and the Eqs. (6.3) in which case the system reduces to only on 
equation: 

(1 - "[,fe) jL,k - (1 - "R.fc) JL,fe + l + "R.fc V?^ = . (6.4) 

Analogous expressions are obtained for the energy injection rates. 

It is here that the non-gradient nature of our models is visible in a nutshell. 
For (6.4) to be of gradient type one needs to have oil^k — '^r k- we have seen 
in Fig. 4, this is, in general, not the case. 

6.1 Infinite volume limit 

We take the continuum limit of Eq. (6.4), with cells of size 1/iV so that the rescaled 
variable ^ — k/(N + 1) is in the domain [0, 1]. Note that the currents cp^ and ip^ 
scale (for fixed external forces) like 1/N. 

In order to close the balance equation for one substitutes the ansatz (6.2) 
into Eq. (6.4). This leads to an involved, but in principle straightforward system of 
nonlinear differential equations, which we do not write down. However, we will 
deal with the special simpler case in which q(Q is approximately constant along 
the chain. 

When q is constant, one finds that T oc 1/j and therefore the assumption (6.2) 
can be reformulated for the corresponding £l and £r to first order in 1/iV as 

J tj JR,k Jl^,k J J J JR,k Jl^,k 

k — ~' ' ' ^R k — ~* * ' (6 5) 

jR,k + JL,fe ' jR,k + 
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Figure 6: 



Dependence of the left (squares) and right (circles) reflection probabilities 
a" — Q!q and on the gradients of the local fields for a chain of 20 RDM cells with 
7 = 0.08. The external gradients were fixed to = 1, = 6, and = = 1500. 
The values for the have been taken from the data of Fig. 3. 



where = —Sq are constants determined by the boundary conditions of (6.4). 
The minus sign in the equation for ^, is a direct consequence of (6.1). 

We have numerically corroborated the hypothesis (6.5). In Fig. 6 we show the 
dependence of the asymmetries and on the gradients of the local injection 
rates for a chain of 20 RDM cells. The asymmetries are seen to depend linearly 
on the local gradient of j. Farther from equilibrium (not shown), deviations will 
appear. 

Using the particle current we express the asymmetry (6.5) in terms of the Jl 
alone as 

J _ -jiL,fc+l 3L,k ^ V^^ 

JL,k+l ^ JL,k V ,^ ^, 

J (6.6) 

J _ ^J-^L,fc+l JL,k 

k — ~ '• J * 

JL,fc+l + JL,k ~ f 

Finally, inserting (6.6) into (6.4) we obtain a closed equation for Jl with only 
one system dependent free parameter, A^. In the scaling limit one finds the fol- 
lowing nonlinear differential equation for the injection rates: 
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Figure 7: Numerical verification of Eq. (6.7), for for systems of size = 20. The 
external parameters are those of Fig. 6. The circles correspond to the r.h.s. of (6.7) as 
explained in the text, while the squares correspond to (jl — Jl k)(^ + 1) — measured 
at ^ = k/{N + 1) — which is an approximation to j{^{0- 



The same equations hold for the "R" versions but with different boundary condi- 
tions. 

Note that Lp^ is the current of a chain of length A^, and hence is asymptotically 
equal to /N. Using (6.7), can be determined by the boundary conditions as 

where W\^) is the coefficient of (p^ in Eq. (6.7). 

As discussed in Sec. 5 (see e.g.. Fig. 3), we have found that at sufficiently 
low or high densities, the a depends little on the local fields, taking practically a 
constant value. Therefore, if the external gradient imposed by the injection of the 
reservoirs into the system is such that the density in all the cells of the chain is 
either low (n < 0.3) or high (n > 10) then the aQ{n{^)) can be taken as constant 
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along the chain. In this situation Eq. (6.7) gives a linear profile for the injection 

rate, S^JlCO = const. In particular, in the high density regime, linear profiles 
will be obtained even when the gradients are very large and the system is far from 
equilibrium. 

For intermediate densities, it follows from (3.4) that ckq is a function of 

/ Ol + Jr)' 

which to lowest order in l/A?^ is 

«(Occy||. (6.8) 

Therefore, in the intermediate density regime, where we have found that aQ(n) 
is a linear function of n, Eq. (6.7) leads to a system of two explicit coupled differ- 
ential equations for j and q, which can be solved. 

The consistency of our approximations was checked for the RDM model by 
combining the measurements of all the quantities appearing in Eq. (6.7). The 
results for the discrete version of this equation are summarized in Fig. 7, which 
shows good agreement with the theory. These measurements also show, that while 
correlations are to be expected for high gradients, [5, 10], they do not seem to 
affect the validity of our approximations. 

The external parameters used are those of Fig. 6, which guarantee that all den- 
sities along the chain lie in the domain of linear dependence shown in Fig. 3, i.e., 
between 0.3 and 10. The data in Fig. 7 were obtained as follows: The derivatives 
Jl were approximated with centered differences along the profiles. The q.q was 
read off the linear fit of Fig. 3. The flux ip is directly read off the simulations, 
while the A coefficients are obtained from the slope in Fig. 6. 

7 Conclusions 

Guided by a stochastic description, we developed a phenomenological description 
of heat and matter transport for a family of Hamiltonian models of many interact- 
ing particles. The description is based on a careful analysis of the way in which 
particles and energy enter and leave the individual cells. 

A stochastic model for a Hamiltonian system must include an appropriate 
treatment of the memory effects that arise due to the deterministic character of 
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the dynamics. This approach would be cumbersome. Instead, we followed a dif- 
ferent approach that, inspired by an approximated stochastic model, is based on 
the phenomenology of reflection probabilities of particles and of energy in terms 
of the local fields and the local gradients. 

We have expressed a set of balance equations, accounting for particle and 
energy conservation in the steady state in terms of the reflection probabilities. 
After specifying a phenomenological law for these probabilities we obtain closed 
expressions for the local fields and show that they capture the essential features of 
the microscopic dynamics including memory effects inherent to any Hamiltonian 
deterministic system. 

A useful observation was to identify two contributions to the memory effects: 
one of geometrical and another of dynamical origin. Close to equilibrium the 
geometric component dominates and the corresponding correction preserves the 
"gradient" character of the system. The theory remarkably predicts the transport 
properties of our class of Hamiltonian models, namely the fluxes of heat and mat- 
ter and the steady state energy and density profiles. 

Far from equilibrium, the dynamical effects strongly depend on the local gradi- 
ents of the thermodynamical fields and thus, on the particular interaction between 
particles and energy tanks. Finally, we have shown that in the continuum limit 
the system is no longer of "gradient" type and the energy and particle currents 
are not proportional to the external gradients. We have obtained the lowest-order 
deviation from the theory in [1]. However, we have found that the gradient type 
condition is restored in the limit of very high or very low densities. 

Given that the dynamical memory effects depend on the particular nature of 
the interaction, our explicit solution (e.g., the particular case of (6.7)), cannot be 
applied to a general context. However, the program outlined and the variables 
used in our derivation are valid and appropriate for the family of models that we 
have considered. 
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